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Abstract 



A model of a cloud formed by massive strings is used as a source of Bianchi type II. We assumed that the 
expansion (9) in the model is proportional to the shear (a). To get exact solution, we have considered the 
equation of state of the fluid to be in the stiff form. It is found that the bulk viscosity plays a very important 
rule in the history of the universe. In presence of bulk viscosity the particles dominate over strings whereas 
in absence of it, strings dominate over the particles which is not in consistence with the recent observations. 
Also we observe that the viscosity caused the expansion of the universe to be accelerating. Our models are 
evolving from an early decelerating phase to a late time accelerating phase. The physical and geometrical 
' (— I ' behavior of these models are discussed. 
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1 Introduction 

• i— i ' 

In recent years, there has been considerable interest in the study of the role of cosmic strings in the evolution of 
the universe, specially at it's early epoches. One of the predictions of Grand Unified Theories (GUT) is that the 
universe underwent a phase transition as the temperature falls down below the Tqjjx ^ 10 28 i"T when the age of 
universe was t GUT - lCT 36 s (Zel'dovich et al. 1975; Kibble. 1980; Everett. 1981; Vilenkin. 1981). There is a 
loss of symmetry when the universe undergoes the GUT phase transition at tcuT- At T < Tqut, the symmetry 
between the strong and electro-weak forces spontaneously broken. The phase transitions associated with loss of 
symmetry leads to the formation of topological defects such as domain walls, cosmic strings, monopoles, etc. As 
mentioned by Kibble (1976), cosmic strings which are important topologically stable defects, might be created 
during a phase transition in the early universe. Although at the present time the existence of the cosmic strings 
can not be detected through our observations, their existence in the early epoches of cosmic evolutions is a 
well established fact. Zel'dovich (1980) was believed that the vacuum strings may generate density fluctuations 
sufficient to explain the galaxy formation. Since the cosmic strings have coupled stress-energy to the gravita- 
tional field, the study of gravitational effects of such strings seems to be interesting. The general relativistic 
treatment of strings was initiated by Letelier (1979, 1983). Here we have considered gravitational effects, arisen 
from strings by coupling of stress energy of strings to the gravitational field. Letelier (1979) defined the massive 
strings as the geometric strings (massless) with particles attached along its expansions. 

As noted by Bclinchon (2009), the string that form the cloud are the generalization of Takabayasis relativistic 
model of strings (called p-strings). Therefore, a cloud model of strings is a model in which we have particles and 
strings together. However, since there is not any observational evidence for the exitance of the cosmic strings 
at the present time, one can eliminate the strings and end up with a cloud of particles at the present time of 
evolution of the universe (Banerjee et al. 1990; Yadav et al. 2007; Saha and Visinescu. 2008, 2010). 



Most cosmological models assume that the matter in the universe can be describe by 'dust' (a pressure-less 
distribution) or at best a perfect fluid. To have realistic cosmological models we should consider the presence of 
a material distribution other than a perfect fluid. Cosmological models of a fluid with viscosity play a significant 
role in the study of evolution of universe. The viscosity mechanism in cosmology can account for high entropy 
per baryon of the present universe (Weinberg. 1972). It is well known that at an early stage of the universe 
when neutrino decoupling occurs, the matter behaves like a viscous fluid (Kolb and Turner. 1990). Weinberg 
(1971, 1972) derived general formulae for bulk and shear viscosity and used these to evaluate the rate of cosmo- 
logical entropy production. He deduced that the most general form of the energy-momentum tensor allowed by 
rotational and space-inversion invariance, contains a bulk viscosity term proportional to the volume expansion 
of the model. Padmanbhan and Chitre (1987) have also noted that viscosity may be of relevance for the future 
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evolution of the Universe. If the coefficient of bulk viscosity decays sufficiently slowly, i.e., £ ~ p n , n < i, then 
the late epochs of the Universe will be viscosity dominated, and the Universe will enter a final inflationary era 
with steady-state character. Cosmological models with viscous fluid in early universe have been widely discussed 
in the literatures (for example see Pradhan et al. 2012; Pradhan and Lata. 2011; Pradhan. 2009; Pradhan and 
Shyam. 2009; Pradhan et al. 2008; Pradhan et al 2007; Pradhan et al. 2005Yadav. 2011 a, b; Yadav et al. 2012). 

The equation of state of a cosmic fluid is defined as the ratio of it's pressure and energy density Lew = 2. 
A fluid in which p — p, is called "stiff fluid" . In this case, the speed of light is equal to speed of sound and 
its governing equations have the same characteristics as those of gravitational field (Zel'dovich. 1970). The 
relevance of stiff equation of state p = p to the matter content of the universe in the early state of evolution of 
universe has first discussed by Barrow (1986). An exact solution of Einsteins field equation with stiff equation 
of state has investigated by Wesson (1978). Mohanty et al. (1982) have investigated cylindrically symmetric 
Zeldovich fluid distribution in General Relativity. Gotz (1988) obtained a plane symmetric solution of Einsteins 
field equation for stiff perfect fluid distribution. Pradhan and Kumhar (2009) have investigated LRS Bianchi 
type II bulk viscous universe with decaying vacuum energy density in General Relativity. Recently A. K. Yadav 
et al. (2011) have investigated string L.R.S Bianchi type II universe in general relativity. 

Bianchi type II space-time has a fundamental role in constructing cosmological models suitable for describing 
the early stages of evolution of universe. Asseo and Sol (1987) emphasized the importance of Bianchi type II 
universe. In the present paper we have considered a locally rotationally symmetric (LRS) model of spatially 
homogeneous Bianchi type-II cosmology. To obtain exact solutions, the field equations have been solved for the 
case when the equation of state of the fluid is in the stiff form. The paper is organized as follows. The metric 
and the field equations are presented in Section 2. In Section 3, we deal with solution of the field equations with 
cloud of strings. In Subsect. 3.1, we describe some physical and geometric properties of the model . In Subsect. 
3.2, we give the solution in absence of bulk viscosity. A dark energy interpretation of the derived models is given 
in Section 4. Finally, in Section 5, concluding remarks are given. 



2 The Metric and Field Equations 

We consider the Bianchi type II metric in the form 

ds 2 = -dt 2 + B 2 (dy + xdzf + A 2 (dx 2 + dz 2 ) (1) 

where A , B are functions of t only. The energy-momentum tensor for a cloud of strings in presence of bulk 
viscosity is taken as 

T ij = (P + P)UiUj +pgi 3 - XxiXj + £0(gi 3 + UiUj), (2) 
where Ui and Xi satisfy condition 

UiU 1 = —x 1 Xi — —1, u l Xi — 0, (3) 

p is the isotropic pressure, p is the proper energy density for a cloud string with particles attached to them, A 
is the string tension density, u l the four-velocity of the particles, and x % is a unit space-like vector representing 
the direction of string. In a co-moving coordinate system, we have 

= (0,0,0,1), a:* = (1,0,0,0). (4) 

The particle density of the configuration is given by 

p p = R ij u l u 3 , (5) 

where p p is the rest energy density of the particles attached to the strings. The string tension density,A, can 
take positive or negative values. Negative value of A represents a universe filled with no strings but only an 
anisotropic fluid whereas its positive value represents strings loaded with particles forming the surface of world 
sheet [8]. 

The Einstein's field equations (with 8irG — 1 and c = 1) 

Rij ~ \ R 9ij = - T *3 ( 6 ) 
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for the metric ([1} leads to the following system of equations: 

A A 2 3 B 2 

G22=2- + ---— = + + A (7) 
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where an overdot stands for the first and double overdot for second derivative with respect to t. 
The spatial volume for LRS B-II is given by 

V = A 2 B. (10) 
We define S = {A 2 B)s as the average scale factor of LRS B-II model (JTJ so that the Hubble's parameter is given 

by 

TT s i (iA b\ 

H= S = 3{-A + b)- 
We define the generalized mean Hubble's parameter H as 

H =^(H X +H y + H z ), (12) 

where H x = 4, H y = ^ and H z — H x are the directional Hubble's parameters in the directions of x, y and z 
respectively. 

The deceleration parameter q is conventionally defined by 

The scalar expansion 6, shear scalar a 2 and the average anisotropy parameter Am are defined by 



2A B 
~A + B' 



(14) 



A m — 



1—1 



where Ai^ = Hi - H(i = 1, 2, 3). 
The Raychaudhuri equation reads 



3| = -2^ + |^-^(p + 3p). (17) 



3 Solution of the Field Equations 

The field equations (|ll[) - (|13p are a system of three equations with five unknown parameters A, B,p, p, A. two 
additional constraints relating these parameters are required to obtain explicit solutions of the system. Firstly, 
we assume that the expansion 6 in the model is proportional to the shear a. This condition leads to 
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which yields to 



A B 
— = m — 
A B 



where m = 2 ^^^ an d b are constants. Eq. (fT9l) . after integration, reduces to 

A = B m , 

Secondly, we assume that the fluid obeys the stiff fluid equation of state i.e 

p = p. 

Using eqs. (|8|) and (|9]) in eq. (|2T|) we get 

A 



B AB A 2 „„ 
A + B +3 AB + A^-^ = °- 



In view of eq. (|20|). eq. (f2"T|) is taken as 



B 2 Y 
25 + 4m — = ; A - 
B (m + 1) 



(19) 

(20) 
(21) 

(22) 
(23) 



where I have assumed that the coefficient of bulk viscosity is inversely proportional to expansion i.e £9 = x (say) 
= constant. 

Let B = f(B) which implies that B = //', where /' = ^. Hence ([23]) can be written as 



X 



B, 



which on integration gives 



dt 



dB 



^JaB 2 + bB- 



-4m 



Here, a 



— — — and b is a positive integrating constant. Hence the model ([T]) is reduced to 
dB 2 



ds z = 



+ B z (dx + zdyf + B zm (dy z + dz z ) 



(24) 
(25) 

(26) 

(27) 



Here we discuss some physical and kinematic properties of string model (|27j) . 

The pressure (p), the energy density (p), the string tension (A), and the particle density (p p ) for the model (|27l) 
are given by 



aB 2 + bB~ 4m 

After using a suitable transformation of coordinates the model (|26|) reduces to 

dT 2 



da' 



aT 2 + bT~ 4 



+ T\dx + zdy) z + T 2m {dy A + dz 2 ). 



3.1 The Geometric and Physical Significance of Model 



P = P = 



m(m + 2) 



(2m 2 + 3m + 1) 



_ 1 /2m 2 - 3m + 1 



2 V 2m 2 + 3m + 1 



2bT~ 2(1+2r 



1 



1 1 



4 T 2(2m-1) ' 



y2(2m-l) ' 



1 / m — 5m + 1 



Pp = ~o 



2 V 2m 2 + 3m 



T )x + -(- + 2)6T- 2 ™ + ^^ ;5I I rrTT . 



(28) 
(29) 
(30) 



From Eqs. (|28l) and ([3U|) . we observe that the energy density p and the particle density p p are decreasing 
functions of time. This behavior of p and p p is shown in figure 1. Also the energy conditions, p > and p p > 
are satisfied under conditions 
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(31) 
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Figure 1: The plot of energy density p, particle density p p and tension density A Vs. T for m = 2, \ — 0-1 an d 
6 = 1 
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respectively. Also A > under 



T> 



2m 2 + 3m + 1 
m 2 — 5m + 1 



(32) 



(33) 



From Eq. (|29p , it is observed that A is an increasing function of time which is always negative and tends to zero 
at late time. It is pointed out by Letelier (1979) that A may be positive or negative. When A < 0, the string 
phase of the universe disappears i.e. we have an anisotropic fluid of particles. This behavior of tension density 
A is also depicted in figure 1. 



To study the behavior of strings and particles in the universe here we define the following parameter 



Pp_ 

|A| 



■3 ( 2^+3^+1 ) X + m(m + 2)6T- 2 ( 1 + 2 ™) + § 



4 T 2(2„ 



2m 2 — 3m + 1 , 
2rn 2 +3m + l / X ' 



(34) 



T 2(2n 



As mentioned before, since strings are not observed at the present time of evolution of the universe, in principle 
we can eliminate the strings and end up with a cloud of particles. In other word, we can say the particles 
dominate over the strings at the present time of the evolution of the universe. Figure 2 is clearly shows that in 
the universe which is describe by the model (|27|) the strings dominate over the particle at initial time whereas the 
particles dominate over the strings at late time. Also it is worth to mention that from figures 2 and 3 we observe 
that at the initial time, when the universe is in the decelerating phase, the strings dominate over the particles 
(p p < A) whereas when the universe is in the accelerating phase particles dominate over the strings {p p > A). 
This is in agreement with the results obtained in Ref (Weinberg. 1976) and (Belinchon. 2009) and also with 
the astronomical observations which predict that there is no direct evidence of strings in the present-day universe. 



The expressions for the scalar of expansion 6, the average generalized Hubble's parameter H, magnitude of 
shear er 2 , proper volume V, deceleration parameter q and the average anisotropy parameter A m for the model 
([27]) are given by 



3H = (1 + 2m) 



2(2m 2 + 3m+l 



(35) 
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Figure 2: The plot of j& Vs. T for m = 2, x = 0.1 and & = 1 
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Figure 3: The plot of deceleration parameter q Vs. T for to = 2, x = 0.1 and 6=1 
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From (pTT|) we observe that 



1 /2o\ 2 < 1+2m > 1 /2o\ 2 < 1+2m > 
q>0 if m> -- &T < I — ) or m< — - & T > I — 

2 \ a J 2 \ a 
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2 V a / 2 \ a 



(40) 



(41) 



A positive sign of q corresponds to the standard decelerating model whereas the negative sign — 1 < q < 
indicates the inflation. Recent observations show that the deceleration parameter of the universe is in the range 
— 1 < q < and the present day universe is undergoing an accelerated expansion. From figure 3 we observe that 
the model (1271) successfully describes the expansion of our universe from decelerating to accelerating phase. 
Also we note that for 

i 

T= (1 - m )x 

[b(2m + 7)(2m 2 + 3m + 1) 

q = — 1 as in the case of de-Sitter universe. 



In absence of any curvature, matter energy density (f2 m ) and dark energy (J7a) are related by the equation 

n m + n A = i, (42) 



where fl m — jjp- and = j^p. Thus equation (|4"2"1) . reduce to 
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A 
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Using equations (l28l) and (|35|) . in equation (|43|) . the cosmological constant is obtained as 



(43) 
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Figure 4: The plot of cosmological constant A Vs. T for m = 2, x = 0.1 and 6=1 
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(44) 



From (I44|) we observe that A is a decreasing function of time and is always positive for to > —0.5 and m < — 1. 
This behavior of cosmological constant A is clearly depicted in figure 4. Recent cosmological observations suggest 
the existence of a positive cosmological constant A with the magnitude A(Gh/c 3 ) « 10 -123 . These observations 
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on magnitude and red-shift of type la supernova suggest that our universe may be an accelerating one with 
induced cosmological density through the cosmological A-term. Thus, our model is consistent with the results 
of recent observations. 

It is worth to mention that for m = 1 from (|4"4"|) we find 

This supports the views in favor of the dependence A oc T~ 2 first expressed by Bertolami (1986 a, b) and later on 
observed by several authors (Rahman. 1990; Chen and Wu. 1990; Berman. 1990 a, b; Berman and Som. 1990; 
Pradhan and Kumar 2001). A relation like equation (jJSJ) also can be found in Brans-Dicke theories when one 
supposes variable gravitational and cosmological "constant" (Peebles and Ratra. 2003; Carmeli and Kuzmenko. 
2002; Gasperini. 1987). We have derived the same variation of A with time in string viscous cosmology in this 
article. 

From the above results, it can be seen that the spatial volume is zero at T — and it increases with the 
increase of T . This shows that the universe starts evolving with zero volume at T — and expands with cosmic 
time T . In derived model, The energy density p, particle density p p , tension density A and the cosmological 
constant A become zero at T — > oo and tend to infinity at T = 0. The model has the point-type singularity 
at T — (MacCallum. 1971). The expansion scalar and shear scalar all tend to zero as T — > oo. The mean 
anisotropy parameter is uniform throughout whole expansion of the universe when m ^ —\ but for m = — \ 
it tends to infinity. This shows that the universe is expanding with the increase of cosmic time but the rate of 
expansion and shear scalar decrease to zero and tend to isotropic. Since % ^constant provided m ^ — |, the 
model does not approach isotropy at any time. But for m = 1 our solution provide a totally isotropic universe. 

3.2 Solutions in Absence of Bulk Viscous 

In absence of bulk viscosity, i.e. x ^ or a - > 0, the metric (|26p reduces to 

rlT 2 

ds ' = -yf^ + T {dx + zdyf + T 2m (dy 2 + dz 2 ). (46) 



The pressure (p), the energy density (p), the string tension (A), and the particle density (p p ) for the model ((46]) 
are given by 

p = P = m{m + 2)bT- 2 ^ - \j^- y (47) 



1 



y2(2m-l) ' 



(48) 



Pp = m(m + 2)bT- 2 ^ + Jj^^y. (49) 

From eqs. (|4T|) - (|4U1) we observe that p, A and p p are decreasing functions of time and A is always negative. 
The energy conditions, p > and p p > are satisfied under 

T> (46m(m + 2))3 , (50) 

and 

m > km > -2 or m < &m < -2. (51) 

respectively. The behavior of p, A and p p are clearly depicted in Figures 5 as a representative case with ap- 
propriate choice of constants of integration and other physical parameters using reasonably well known situations. 



From eqs. (j48|) and (|49j) we obtain 

y^r = &m(m + 2)T" 4 + 3/4. (52) 
1^1 

Equation ((52|) obviously shows that pq- is a decreasing function of time i.e as time goes on, the strings dominate 
over the particles which in contradict with the result obtained in the first case in presence of bulk viscosity. This 
result is of course is not consistent with the astronomical observations, which predict that there is no direct 
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Figure 5: The plot of energy density p, particle density p p and tension density A Vs. T for m = 2 and 6=1 
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Figure 6: The plot of Vs. T for m = 2 and b = 1 
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evidence of strings in the present-day universe. Therefore, we conclude that the bulk viscosity may plays an 
important rule in creation of the particles from strings. 



The expressions for the scalar of expansion 6, the average generalized Hubble's parameter H, magnitude of 
shear er 2 , proper volume V, deceleration parameter q and the average anisotropy parameter A m for the model 
(|lrJ|) are given by 

9 = 3H = (1 + 2m)VbT-( 1+2m \ (53) 

ff2 = (m- l) 2 bT _ 2(1+2m); (54) 
3 

V = T 2m+ \ (55) 

■ 1 \ 2 
1 — m 



A m = 2 — — - . (57) 
\ 1 + 2m / 

From (|56|) we observe that g > if m > — | and g < if to < — |. But from eq. (|55t we observe that, m < —\ 
represents an accelerating collapsing universe with high blue shift. Since the recent observations (Riess et al. 
2001) indicate that we live in an accelerating expanding universe with red shift, we conclude that in absence of 
bulk viscosity a universe with decreasing rate of expansion is the only possible scenario. 



4 Dark Energy Interpretation of The Models 

Fig. 3 clearly shows that the presence of bulk viscosity in the cosmic fluid causes a decelerating to accelerating 
expansion of the universe. Also from Raychaudhuri's equation. (|17p. we observe that that bulk viscosity can 
play the role of an agent that derive the present acceleration of the Universe. 
In Eckart's theory (1940) a viscous pressure is specified by 



peff — p 



n. 



(58) 



Here II = — is the viscous pressure. Therefore in our models the effective pressure (stiff fluid plus viscous 
fluid) can be written as 



e ft to(to + 2) 

P n =p-x = — z — 



X 



(2m 2 + 3to+ 1) 



26r -2(l+2r, 



1 1 



4 J>2(2m-1) 



(59) 



Using eqs. (|28l) and (I59| the effective equation of state of the net fluid is obtained as 

X 



p 



m(m+2) 
2 



X 



(2m 2 +3m+l) 



2bT-2(l+2r 



4 T 2(2m-1) 



(60) 



The behavior of effective equation of state, uj e ^ , in terms of cosmic time T is shown in Fig. 7. It is observed 
that the w 6 -' * parameter is a decreasing function of T and the rapidity of its decrease depends on the value of %. 
We see that in absence of bulk viscosity the models do not exhibit accelerating expansion (solid line) whereas 
in presence of viscosity our models exhibit a decelerating to an accelerating expansion. From both eq. (|60[) and 
Fig. 7 we observe that at the later stage of evolution, effective equation of state tends to the same constant 
value i.e w e// = 1 - 2(2 ^+^ ) +1) independent of the value of \- 



5 Concluding Remarks 

In this paper we have presented a new exact solution of Einsteins field equations for LRS Bianchi type II space- 
time with cloud of strings which is different from the other authors solution. In general the models are expanding, 
shearing and non-rotating. It is found that in present of bulk viscosity particles dominate over the strings at late 
time whereas in absence of viscosity strings dominate over the particles which is a contradictory result. In other 
hand, from Raychaudhuri's equation (|17|) we observe that that bulk viscosity can play the role of an agent that 
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Figure 7: The plot of effective equation of state uj e ^ Vs. T for m = 1 and 6=1 

derive the present acceleration of the Universe. Hence we conclude that the bulk viscosity plays an important role 
on the evolution of the universe. For a universe which was decelerating in past and accelerating at the present 
time, the deceleration parameter must show signature flipping (see the Refs. Padmanabhan and Roychowdhury. 
2003; Amcndola. 2003; Caldwell et al. 2006). Our models are evolving from an early decelerating phase to a 
late time accelerating phase (see, Figure 3) which is in good agreement with recent observations (Riess et al. 
2001). 
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